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We propose a mechanism for the fast dissipation of magnetic fields which is effective in a stratified 
medium where ion motions can be neglected. In such a medium, the field is frozen into the electrons 
and Hall currents prevail. Although Hall currents conserve magnetic energy, in the presence of 
density gradients, they are able to create current sheets which can be the sites for efficient dissipation 
of magnetic fields. We recover the frequency, ujmh, for Hall oscillations modified by the presence 
of density gradients. We show that these oscillations can lead to the exchange of energy between 
different components of the field. We calculate the time evolution and show that magnetic fields can 
dissipate on a timescale of order 1/umh- This mechanism can play an important role for magnetic 
dissipation in systems with very steep density gradients where the ions are static such as those found 
in the solid crust of neutron stars. 
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I. INTRODUCTION 

Rapid dissipation of magnetic fields is currently one of 
the key problems in astrophysics. On account of the gen- 
erally large electrical conductivities that obtain in astro- 
physical settings, the Ohmic dissipation of fields usually 
takes place on very long time-scales. However, it is quite 
often observed that astrophysical magnetic fields change 
topology on a very short time-scale, giving rise to a vari- 
ety of transient phenomena. An explanation of such fast 
changes is crucial to the understanding of solar activity, 
in particular, solar flares H], and other active phenom- 
ena observed in stars [y||]T It is also known that fast 
reconnection of magnetic fields is basic to the operation 
of non- linear dynamos ^ . 

If a current sheet is formed in a plasma, the recon- 
nection takes place slowly due to the time-scale for the 
removal of matter from the site of reconnection, as in 
the Parker-Sweet mechanism Indeed, if the dissipa- 
tion in the current sheets were to be fast, with the field 
moving with Alfven speed toward the current sheet and 
getting dissipated there due to reconnection, then the 
matter that is frozen into the field would also move with 
the same speed towards the sheet. The plasma would, 
therefore, accumulate at the sheets and halt the recon- 
nection, unless there is some efficient evacuation process 
operating at the sheet. The sheets are usually narrow 
and the outflow is rather inefficient, even if it takes place 
at the Alfven speed 

We propose a mechanism of fast dissipation of mag- 
netic fields that occurs at modified Hall frequencies. The 
mechanism is relevant in all situation when Hall currents 
predominate, and there is a density stratification. In this 
case, the magnetic field follows the (electric) drift veloc- 
ity of the electrons. In the presence of density gradients, 
the profile of magnetic field changes in such a way that it 
forms a current sheet. This steepening of the front is not 



accompanied by the flow of plasma towards the sheet, 
the drift velocity being parallel to it. Consequently, the 
current sheet is formed and the field is efficiently dissi- 
pated with no accumulation of material in contrast to the 
Parker-Sweet reconnection mechanism. 

The modified Hall frequency that we recover below, 
^MH, occurs in a stratified medium when the ions re- 
main static. Two example of such situations are: pene- 
tration waves in low collision plasmas relevant for plasma 
switches and neutron star crusts. In the case of pene- 
tration waves, the ion response time is long compared to 
the wave timescale and the ions are approximately static. 
In the case of neutron star crusts, the ions form a very 
high density lattice of iron rich nuclei with densities vary- 
ing from ~ 10®g/cm~'^ to ~ I0^^g/cm~'^ under 0.8 km 
§. In both cases, the field dynamic is governed by the 
electron drift motion. 

In the following sections, we discuss how Hall cur- 
rents in a stratified medium can generate fast dissipa- 
tion in the non-linear regime. In Section II, we describe 
the linear Hall oscillations and discuss how poloidal and 
toroidal fields exchange energy during these oscillations. 
We also recover the modified Hall frequency for a strat- 
ified medium. We discuss the limitations of the oscilla- 
tory solutions about a stationary configuration in Sec. 
III. We argue that in general there is no stationary con- 
figuration for large scale magnetic fields and that cur- 
rent sheets develop. The oscillations can occur only "lo- 
cally", i.e., on small scales. In Sec. IV, we show that 
the magnetic field evolution is governed by a non-linear 
equation similar to Burgers equation. We solve the evo- 
lution for a toroidal field configuration numerically and 
show that current sheets develop and magnetic dissipa- 
tion is efficient. The dissipation timescale is ^ 1/lumh- 
We describe numerical solutions for two configurations: 
a toroidal magnetic field of one polarity; and a toroidal 
magnetic field consisting of two oppositely directed fields. 
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We show that these fields evolve towards forming current 
sheets that rapidly dissipate. In Sec. V, we relate the 
dynamics of toroidal fields with that of poloidal fields, 
and summarize the different possibilities of the field evo- 
lution. We close by discussing the application of this 
physical mechanism focusing particularly on the case of 
neutron stars crusts (Sec. VI). 



II. LINEAR OSCILLATIONS. 

Consider the magnetic field evolution in the case where 
the motion of ions can be neglected. This is the case for 
neutron stars' solid crust. We will consider coUisional 
plasma. Then, the field evolution follows from Ohm's 
law. 



47r 1 
V X B = — cr(E + -Ve X B), 

c c 



(1) 



where Ve is electron velocity. As the conductivity is usu- 
ally high, the left-hand side of (0) can be neglected, re- 
sulting in 



E 4- -Ve X B = 0, 
c 



(2) 



corresponding to the electric drift of electrons. 

Taking Vx of equation (|l]), we recover the induction 
equation, with the Hall effect. 



V X [ve X B] - V X ?7V X B 
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where rj = (? j^-na. Equation (|^) results in the following 
energy balance, 



2dt 



J B^dV ^- J r]{\/ X B)^dV. 



(4) 



In order to describe the evolution of the field in a strat- 
ified medium, we consider and axisymmetric magnetic 
field, which can be expressed as the sum of poloidal and 
toroidal components. 



B — B,] 



In order to simplify the geometry we assume that the ra- 
dius of the star, R, is large compared to the wavelengths 
involved such that we can work in Cartesian coordinates 
on the surface of a sphere. We define the x and y axes 
in the horizontal plane as the latitudinal and azimuthal 
(longitudinal) directions respectively, while z is the ver- 
tical direction. Then, the poloidal field is described by 



Bp = {Bx{x, z), 0, B:,{x, z)}, 



while the toroidal field is given by 

Bt = {0,By{x,z),0}. 

If the resistivity rj can be neglected, which is justified 
in highly conducting media such as neutron star crusts, 
then the magnetic energy is conserved, according to (Q). 
Therefore, all that happens to the magnetic field are os- 
cillations about a stationary configuration. 

Consider, for example, an initial poloidal field. Bp = 
{So, 0,0}, where Bq is a constant background field. As- 
suming first that the density is also constant, and con- 
sidering small perturbations of the magnetic field of the 
form, 



(5) 



{ky — because of axial symmetry), we find substituting 
in H for r/ ^ 



|(k-a;e)|c2fc 



(6) 



where We = gBq / mc is the electron cyclotron frequency, 
and LOp is the plasma frequency. We have thus recov- 
ered the well known Hall oscillations or whistlers. Note 
that even if &j, = initially, i.e., the toroidal compo- 
nent is absent, it will be generated reaching the level of 
the (perturbed) poloidal component; and thus, the en- 
ergy will be exchanged between the poloidal and toroidal 
components. 

The situation is different if the large scale background 
field is toroidal, i.e., Bq = {Q,Bq,Q}. Then a perturba- 
tion of the form (H) would not result in oscillations (|^) , 
because (k • We) — 0. 

Let us now recall that the density is not a constant, 
but, rather, it has a steep dependence on z. By including 
the spatial dependence of the density in (^ , we recover 
the Hall frequency modified by the presence of density 
gradients. 



^MH — 



(k • [we X Vn])c^ 



(7) 



Note that is a pseudo- vector, and therefore the fre- 
quency ujmh is a real scalar, as it should be, just as loh 
is a real scalar as well, see (^. The phase velocity, cor- 
responding to (^ can be written as 



^MH — 7 



X We 



(8) 



and ujmh = k • "Vmh ■ A similar case is known in low col- 
lisional plasmas where the corresponding wave is called 
magnetic penetration wave ||^. 

The wave described by (0,||) corresponds to only 
toroidal perturbations due to the chosen initial configu- 
rations. In this special case there is no poloidal field ini- 
tially and no energy exchange occur between toroidal and 
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poloidal components. However, in general the two com- 
ponents are present and this exchange does take place. In 
order to see this, let us return to the large scale poloidal 
field, Bo = {Bo,0,0}, taking into account that the den- 
sity is a function of z. We look for solutions of the lin- 
earized equations in the form. 



b = {dza{z), by{z), -ikxa{z)}e 



(9) 



cf. (|^). Then we obtain the following dispersion relation, 

(10) 



-g {k^-d,d,)a. 



In order to estimate the frequency consider two zones 
Z2 < z < 2i, with density n2, and zi < z < = 
is the top of the crust), with density rii, and |zi| = /ii, 
and zi — Z2 = /i2, /ii,2 being the scale hight in these two 
zones. Assuming that n2 S> jt-i, and /12 ^ ^i, we obtain 



TT |(k • We) 



-|fcx + l/hi 



(11) 



"^2 



where 



is the plasma frequency based on the den- 



sity 712 ■ Note that if fca, <C (large horizontal length 
scale), the frequency is essentially the same as ujmh in 
(0). This is the main characteristic frequency of mag- 
netic fluctuations in the crust due to the steep density 
gradient. As seen from (^), this mode does involve both 
poloidal and toroidal components. 

The most general case involves non-linear coupling be- 
tween the poloidal and toroidal fields. Qualitatively the 
same situation will take place: if we start with a poloidal 
field, supported by the currents in the crust, a toroidal 
field will be generated. The current velocity is toroidal, 
and, according to (|^), the toroidal field is stretched out 
from the poloidal, analogously to the effect of differen- 
tial rotation. However, unlike the latter, the Hall current 
conserves the energy, and therefore the new toroidal field 
will grow at the expense of the poloidal field. In other 
words, while the strength of the toroidal field is increas- 
ing, that of the poloidal component should decrease. Of 
course, the toroidal field cannot grow indefinitely under 
these circumstances, and eventually the field will either 
reach some steady state, or the poloidal and toroidal 
fields will exchange their energies, oscillating with fre- 
quency UMH- 

Note, however, that including dissipation may drasti- 
cally change the situation, and, in some cases, discussed 
below in Sec. Ill, and IV, the field will rapidly dissipate 
instead of oscillate. 



III. THE PROBLEM OF STATIONARY STATES. 

This simple picture of oscillations implicitly assumes 
that they proceed about some stationary state, which 



presumably exists. The large scale background field con- 
sidered above was uniform and trivially stationary. We 
will show that, in general, the large scale field is not sta- 
tionary but evolves with time. 

It is clear from (||) that the stationary state is possible 
if, neglecting diffusion. 
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[V X B] X B = V$, 



(12) 



that is, the electric field is potential. We will show that 
condition ( p^ ) does not trivially occur even for extremely 
simple topologies, due to the gradient of the density. In- 
deed, consider an initial configuration consisting only of 
a toroidal field. 

Equation (|^) for a pure toroidal field can be written as 

VzdzBy = viV'^By, 



dfK, 



r, dx B), 



where 



Vx = 



dxVi 



(13) 



dzV- (14) 



If we neglect the resistivity in this expression, we recover 
the penetration wave velocity (Bh, v — > v^ih as 77 0. 

It can be seen from equation ( 14) that the a;-component 
of the velocity is non-vanishing, due to the vertical gra- 
dient of the density. Note that both the density gradient 
and the gradient of the resistivity 77, are negligible in the 
X direction, and the Vz component defined only by the 
resistivity gradient is also small. Since any toroidal field 
should vanish at least at the two poles, there is always a 
latitudinal dependence of the toroidal field, that is to say 
that By is always a function of x. Hence, according to 
(|l3|), the toroidal magnetic field can never attain a sta- 
tionary state. In other words, the electric field cannot be 
irrotational, as in (|l^), and its non-potential part results 
in the time evolution of the magnetic field. 

Note that in infinite space equation (|l^) conserves 
magnetic flux. 



Bydxdz — const, 



(15) 



but the magnetic energy is dissipated according to, 

~jBldxdz = -j7]{\7By)^dxdz, (16) 

which is a particular case of (Q). 

On the other hand, for a real toroidal field which should 
vanish at the poles, i.e., at x = ±7ri?/2, the magnetic fiux 
is not conserved. Indeed, according to (|l3|), 

J Bydxdz ^ J ridxBy{x ^ TTR/2)dz 



- J 7]d^By{x = -■nR/2)dz 



(17) 



The right hand side gives considerable contribution when 
current sheets are formed at x = ±.ttR/2. 
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IV. TOROIDAL MAGNETIC FIELD 
EVOLUTION: FORMATION OF CURRENT 
SHEETS. 

A. Analytical and Numerical Solutions. 

In order to study the evolution of the field, according to 
( p^ ) , we reduce this equation to (neglecting the resistivity 
gradient, and resistive diffusion in z-direction) 



dtb + hdxb = rjdxdxh, 



where 



b = Byp, 



P ■ 



A'Ke'n? 



(18) 



(19) 



This is, in fact, the Burgers equation, the exact solution 
of which is well known, see, e.g., First, let us illustrate 
a solution in the form of a traveling shock wave. 



1 — tanh 



{x - bat)bo 
2r, 



(20) 



where bo is a constant, cf. e.g. [|6|j9|]. The penetration 
wave ( pO| ) does not decay because the magnetic field is 
pumped into the system from — oo. Therefore, it is more 
appropriate for our purposes to use the general exact 
solution, which we recover by using the transformation 



b = -2ridx In^, 



(21) 



to get. 



(47r?7i)i/2 



g-(.-.')V(4r,t)^(^^0,:E',z)dx' . (22) 



Generally, the toroidal field By is a function of both x and 
z, and, since the z-dependence enters only parametrically 
into (|l^ , [l9| ), the solution (|2l|j2^ ) can be used for each 
level z — const. 

To illustrate the time evolution of the magnetic field, 
we demonstrate the following two cases. In the simplest 
case, we assume that the toroidal magnetic field does 
not change sign. Then the horizontal velocity Vx is ex- 
pected to drive the field to one of the poles, either to 
the South, or to the North, depending on the sign of the 
field. The gradient of the field steepens, as in a shock 
wave, thus forming a current sheet, where the magnetic 
field is finally dissipated. In the second case, consider two 
toroidal fields with opposite polarities in the two hemi- 
spheres. The toroidal field vanishes at the equator. We 
then expect that the two toroidal fields can be driven 
by the latitudinal velocity towards the equator, where 
the current sheet is formed, and the fields are efhciently 
destroyed. 



The integral in (22) was calculated numerically, and 
then the distribution of magnetic field By was recovered 
from equations ([l^) and (|2^). Let us discuss the first 



case where the toroidal magnetic field does not change 
sign. Its evolution is depicted in Fig. 1, where the initial 
field distribution is indicated by the dashed line. The 
field profile starts to steepen in very few turn-over time 
steps, and moves towards the polar region. Note that the 
magnetic field in Fig. 1 is not pumped into the system, 
cf. (pO[). Therefore, unlike the traveling wave (pO|), as 
the magnetic field spreads its amplitude decreases keep- 
ing the magnetic flux conserved and, thus, the same area 
under each curve, see (|l^). As a result of decreasing 
magnetic field, the process slows down, because the pen- 
etration velocity (^) is proportional to B, and therefore 
it decreases as well. 




FIG. 1. Evolution of magnetic field of single polarity in 
the crustal region of a neutron star, as it approaches a polar 
region. The equator is located at zero latitude, and the time 
t is expressed in units of turn-over time to, given in (p^). 

In infinite space, both the shock wave ( pO| ) and the so- 
lution depicted in Fig. 1 do not result in dissipation of 
magnetic field and the magnetic flux is conserved accord- 
ing to (p^. The field is only spread out. However, for a 
finite case such as that of a star, the boundary conditions 
at the poles forces the field to go to zero. When the shock 
wave reaches the pole, a current sheet is formed and the 
field starts to dissipate according to (|l^) . Eventually, the 
magnetic flux goes to zero. 

In order to see this dissipation at a zero-point, we 
proceed to the second example. Namely, consider the 
toroidal field changing sign at x = 0. It is straightfor- 
ward to construct a solution, analogous to the traveling 
wave (EG), 



b = — 6n tanh 



xbo 
2ri 



(23) 



Similarly to the Parker-Sweet solution Q, the magnetic 
field of opposite polarities is transported from x — > ±oo 
with "velocity" ±6o, and is dissipated at a; = 0. The 
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solution is stationary because the boundary conditions 
are: By{x ±00) ±bo/p, see (19). 




These analogies, although helpful, cannot be taken 
completely due to the topological constraints on mag- 
netic fields. For instance, since the magnetic structures 
are frozen into the electron fluid, they are generally more 
persistent than vortices in hydrodynamics. The topology 
of magnetic field cannot be easily changed, a situation 
similarto what one obtains in MHD |^,|2|. The mag- 
netic structures we considered are non-stationary due to 
global effects (like the boundary conditions at the poles) 
and not locally unstable like the case of vorticity and 
magnetic turbulence. Besides, in the case of magnetic 
turbulence, the characteristic frequency coincides with 
the Hall oscillations frequency (^, which is smaller than 
^MH from d^) for situations with sharp density gradients. 
Therefore, our mechanism is more efficient. 



FIG. 2. The evolution of magnetic fields of opposite polar- 
ities in the two hemispheres. In this case, the two fields are 
approaching each other to form a current sheet at the equator, 
where the magnetic energy is efficiently dissipated. 

In Fig. 2 we display the evolution of a toroidal 
field with opposite polarities in each hemispheres and 
By{x = 0) = 0. The field also vanishes at the poles, 
By{x = ±ttR/2) = 0, which makes the solution evolve in 
time (again, in contrast to the infinite space case). Here, 
the two fields get compressed into each other and form 
a sharp gradient of magnetic field at the equator. Note 
that the total magnetic flux is zero and, of course, triv- 
ially conserved. As to the magnetic energy, it decreases 
dramatically because of the very efficient Ohmic dissi- 
pation at the equatorial region. In this region, a current 
sheet is formed in practically only one turn-over time and 
the field dissipates in the same time scale. An analyti- 
cal estimate of the dissipation time is given in the next 
subsection; it illustrates why the dissipation observed nu- 
merically is so efficient. 

B. Physical Interpretation of the Mechanism. 

In order to develop a physical interpretation of the so- 
lutions above, we draw a few analogies. The equation for 
the magnetic field (||) resembles the vorticity equation for 
incompressible hydrodynamics with high Reynolds num- 
bers. Therefore, the modified Hall drift should lead to 
a situation analogous to a magnetic turbulent state . 
Another interpretation of our solutions is the nonlinear 
interaction of different wavenumber Hall oscillations re- 
sulting in an energy cascade to small scales . As a re- 
sult, the field gradients steepen and we get an enhanced 
local rate of Ohmic dissipation which provides an effec- 
tive mechanism for the dissipation of magnetic energy. 




FIG. 3. Origin of the penetration velocity. The magnetic 
field evolution is depicted on the top of the figure. The di- 
rection of the field is out of the page. It can be seen that the 
field profile steepens during the evolution due to the change 
in drift velocity (double arrows). The drift proceeds in the 
vertical direction, and, therefore, there is no accumulation of 
matter at the current sheet. 

In magnetohydrodynamics, fast reconnection encoun- 
ters difficulties because the rapid transport of mag- 
netic fields towards the current sheet, where the energy 
is dissipated, is accompanied by a plasma movement in 
the same direction. The evacuation of matter from the 
current sheet limits the rate of reconnection: as matter 
accumulates and the pressure increases eventually halting 
the movement of the magnetic field towards the current 
sheet, thus preventing further reconnection. The speed 



5 



of the evacuation is limited by the Alfven velocity in a 
narrow current sheet. In our mechanism, we do not en- 
counter this difficulty, because does not transport the 
mass. Indeed, according to ( [l^jl9| ), there is no outflow 
from the current sheet: the magnetic field is transported 
only to the sheet by the penetration velocity v^- This is 
evident from the exact solution ([2^). 

In order to understand why the modified Hall drift or 
penetration velocity does not transport any mass, we first 
note that, generally, the penetration velocity (^) is not 
parallel to the electric drift velocity (||) . So the question 
arises, why the magnetic field is moving in the direction 
of the penetration velocity in the first place? The answer 
is illustrated in Fig. 3. For simplicity, we choose to 
illustrate a field that depends only on the x-coordinate 
(as in (p3|)). The electric drift of electrons (Q) can be 
easily found from Amper's law, 



-VxB, 



(24) 



and, clearly, it proceeds in the vertical direction (depicted 
by double arrows in Fig. 3). It follows from (^J) that 



V • nVe 



0. 



(25) 



Due to the density gradient, the plasma gets compressed 
as it moves down and the descending motion decelerates, 
as follows from (|2^). As a result, the magnetic field am- 
plitude increases. On the other hand, the ascending mo- 
tion is accompanied by the decompression of plasma, and 
correspondingly, the field amplitude decreases. As a re- 
sult, the field profile steepens, as if there were a motion of 
plasma toward the current sheet, that is in the horizontal 
direction. However, as mentioned, the real drift motion 
proceeds parallel to the sheet (in the vertical direction), 
and therefore there is no accumulation of matter in the 
sheet. These circumstances make the fast dissipation of 
magnetic field possible. 

As mentioned above, the penetration wave (0,§) is 
known to propagate as a shock wave in low coUisional 
plasm a [pP ]. Due to the conservation of the magnetic 
flux, (|15|), either the magnetic field is only transported 
by the shock wave (EOh or it just disperses to infinity. 



Only the presence of zero-points result in the destruction 
of magnetic flux. This happens either at the poles or be- 
tween two shock waves with opposite magnetic fields that 
collide. The collisional front width of the shock wave co- 
incides with current sheet thickness. Indeed, the balance 
between the convective and resistive terms in equation 
( p!3| ) appears in a current sheet of thickness. 



Vx 



(26) 



coinciding with characteristic length of both ( |20| ) and 
( p3|) . Recall that the penetration velocity Vx is defined 
in (|14D, and, in a more specific way, in (H). 



It is known that the time scale of magnetic dissipa- 
tion, to, is entirely defined by the velocity with which 
the magnetic field is moving to the current sheet As 
seen from the exact stationary solution (p3|), this speed 
is in effect the penetration velocity Vx- Therefore, 



to = L/vx, 



(27) 



L being the macroscopic latitudinal scale. It is useful 
to confirm this dissipation rate estimate analyzing the 
energy dissipation directly from (16). Namely, we esti- 
mate VBy in ( p^ ) as By/ 5, and the area occupied by the 
current sheet is Sg — SL. Then, we get from (p^). 



to 



q2 



(28) 



from which to is recovered as in (27). Note that the dis- 
sipation time ( |27| ) is independent of the resistivity 77, and 
therefore the process considered here is fast. This ex- 
plains why the dissipation is so efficient in the numerical 
results shown in Fig. 2. 

We finally note that the efficient dissipation depicted 
in Fig. 2 proceeds when the penetration velocities of the 
two toroidal fields point to each other, and therefore they 
collide. If we change sign of magnetic fields, then, accord- 
ing to (H), the penetration velocity changes its direction, 
and, as a result, the two toroidal fields would not collide, 
but instead drift to the polar regions, where they will 
eventually decay. As mentioned above in Sec. IVA, the 
latter process is much less efficient because the magnetic 
field strength decreases as the fields move to the poles, as 
seen from Fig. 1, and therefore the penetration velocity 
decreases. 



V. EVOLUTION OF THE POLOIDAL FIELD. 

In Figs. 1 and 2, we have shown the solution of 
our numerical calculationsfor the evolution of toroidal 
fields with different initial profiles. Since the toroidal 
and poloidal fields are coupled through non-linear oscilla- 
tions, we expect the dissipation of toroidal fields to cause 
the eventual decay of the poloidal field. The exact evo- 
lution of the poloidal field is a harder problem to solve 
at this stage and we leave it for future studies. Below, 
we only discuss the expected qualitative behavior of the 
poloidal field. 

As we saw in Sec. II, the linear oscillations exchange 
energy between the poloidal and toroidal components. 
In other words, an initial poloidal field would generate a 
toroidal one. The generated toroidal field could have the 
nonlinear evolution depicted on either Fig. 1 or Fig. 2. 
In the case that the generated toroidal field is in the con- 
figuration of Fig. 1, it would slowly drift to the poles. We 
expect that, in this case, one would observe oscillations, 
because the dissipation is inefficient. 
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Consider now the case when the toroidal fields are gen- 
erated with the configuration as in Fig. 2. In order to 
follow this generation in the nonlinear case, we write, 
according to (||), 



where 



_9 

dz 



(30) 



is the azimuthal current. Note that, for neutron star 
crusts, the density gradient in the radial direction is ex- 
tremely steep in the crustal regions, spanning some nine 
orders of magnitude over a distance of a few hundred me- 
ters below the surface. The variations of other quantities 
in equation (|2^) can therefore be neglected to get 



d c cdzfi 



dz Anne A-kch? 



jyBz 



(31) 



It is evident from (|31]) that the toroidal field can always 
be generated due to the sharp density gradient, unless 
the poloidal field in the crustal region is current free. If, 
indeed, the field is anchored in the core (meaning that the 
currents supporting the field are confined there), then 
there is no Hall current present. At present, the locus 
where the field is anchored in neutron stars is a matter 
of debate with no clear resolution (see, e.g., ||^,[l3||). 

If we assume that part of the current supporting the 
field is present in the crustal layers, then the toroidal 
field will be generated from the poloidal field by the pro- 
cess outlined above in (|3l|). On the other hand, the to- 
tal magnetic energy is essentially conserved (apart from 
weak Joule dissipation), which means that the newly gen- 
erated toroidal field would result in a back reaction on 
the poloidal field in such a way that the energy of the 
latter is decreased. In the absence of Ohmic dissipation, 
the toroidal field would grow to a certain level, and then 
start to decrease, thus presenting an oscillatory behavior, 
as described earlier at the end of Sec. II. If we incorpo- 
rate Ohmic diffusion in the nonlinear case the both the 
toroidal and the poloidal fields should eventually decay. 

Again, in the case of Fig. 1, the toroidal fields would 
slowly drift to the poles and we expect to observe oscilla- 
tions, because the dissipation is inefficient. In the case of 
Fig. 2, the toroidal field is efficiently dissipated and con- 
sequently, according to equation (^, both the poloidal 
and toroidal fields decay. Indeed, the Ohmic dissipation 
is now increased due to the presence of current sheets, so 
that the equation M) can be written in the form 



ld_ 
2dt 



B^dV 



(32) 



In order to follow the evolution of the poloidal field, we 
introduce efficient dissipation discussed above into equa- 
tion (|3l|), to get 



d By 

TrBy = LUMuBp — — , 
at to 



(33) 



and, because the Hall current conserves the total energy, 
/(Bp -|- By)dV, the back reaction of the toroidal mag- 
netic field on the poloidal component can be expressed 
analogously. 



d_ 

di 



Bp — —LUMHBy. 



(34) 



Seeking solutions ~ e'''*, we find a dispersion relation, 




(35) 



It can be seen from equation (|3g) that the decay time for 
the poloidal component is also of the order of io • 

To summarize, we can delineate three regimes for mag- 
netic fields in the crusts of neutron stars: 

(1) The currents supporting the fields in the crust are 
anchored in the core, i.e., no currents in the crust. Then, 
there is no Hall current (by definition), and no evolution 
of the fields related to the processes we described here. 

(2) The currents or part of the current are situated 
in the crust. Then, the poloidal field inevitably gener- 
ates toroidal fields, which, depending on the sign of the 
initial poloidal field, may result in either penetration ve- 
locity pushing these toroidal fields apart or pushing them 
together. In the first case, the dissipation is less efficient, 
being limited to slow decay at the poles (Fig. 1). Al- 
though slower than the equatorial case, the decay at the 
poles is still faster than the general Ohmic decay. 

(3) In the second case, when the toroidal fields are 
pushed together as in Fig. 2, we expect both toroidal 
and poloidal fields to decay according to (p5[). 

Note that in any of the above cases there would be os- 
cillations on scales small compared with the radius of the 
star (i.e., in the geometric optics limit), with frequency 
^MH ■ These oscillations in general will decay with Ohmic 
decay time. 



VI. DISCUSSION 

Magnetic fields are an important feature of neutron 
stars since, together with the rapid rotation of the star, 
they determine the characteristics of the pulsar emission. 
The source of a wide range of magnetic field strengths 
(~ 10* — 10^^ G) associated with neutron stars is yet to 
be well understood. The seven orders of magnitude span 
may be attributed to the different environments in which 
neutron stars are present: from isolated objects to accrit- 
ing members of binary systems. This range could also be 
the result of different conditions at the time of birth of 
neutron stars, such as the gravitational collapse of the 
progenitor massive star or the accretion-induced collapse 
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of a white dwarf [|13|. In any case, the very high electri- 
cal conductivity renders the Ohmic decay inefficient with 
typical accreting of the order of billions of years. If neu- 
tron star fields decay over their observable lifetime, an 
alternative decay mechanism is necessary to explain this 
behavior. 

One of the uncertainties concerning the evolution of 
neutron star magnetic fields is their location in the stellar 
interior. Should the field be a fossil remnant left over 
from the progenitor star, it could permeate the whole 
body of the neutron star. On the other hand, if the 
magnetic field is generated after a neutron star is born 
via a battery effect or a dynamo process it is likely 
to be confined to its crustal layers. As we discussed in the 
previous section, the exact location of the currents will 
determine if the mechanism proposed here is operating 
in neutron stars or not. 

For instance, the interaction between differential rota- 
tion and magnetic field during the first few seconds of a 
nascent neutron star's life would generate strong toroidal 
magnetic fields in the subsurface layers of the star. With 
the rapid cooling of the star, the crust solidifies with the 
ions forming a lattice in the presence of relativistic elec- 
trons. Some fraction of the toroidal field will have dif- 
ferent signs in the Northern and Southern hemispheres, 
like the one illustrated in Fig. 2. Under these conditions 
the magnetic field is frozen into the electron gas and Hall 
currents in the crustal layers can arise and our mecha- 
nism will be effective. In contrast, the Ohmic dissipation 
in the crustal layers takes place on a very long time scale. 

If part of the currents supporting the fields are situ- 
ated in the crust, we can use our mechanism to estimate 
the timescale for rapid dissipation to occur. Taking typ- 
ical numbers for the crustal layers of a neutron star, at 
density scale height h of lO^cm, n = lO'^^cm"'^ and the 
magnetic field is lO^^G, then ~ 10~^cm/ sec. The cor- 
responding time scale to = L/vx^ where L is horizontal 
scale of the magnetic field, is to — lO^^sec « 3 million 
years, assuming L — lO^cm. On the other hand, for a 
density scale height of 3 • 10'^ cm, then n — lO'^^cm"'^, 
we have for B — lO^^G,?;^; ~ 10~^cm/sec, therefore, 
to = IQ^^sec « 3000 years. Finally, if we take a scale 
height of 10'^ cm, then the density is n = lO"^", and 
Vx ~ 10~^cm/ sec, and to = 30 years. In a real neu- 
tron star, all of these time scales are present if currents 
occur throughout the crust. 

Indeed, due to the sharp gradient of electrical conduc- 
tivity in the crustal region, we can consider the depth 
10'^ cm as a boundary between two layers, with different 
CTi.2, index 1 corresponding to the upper layer, and index 
2 to the lower layer, with a2 3> cri. Then, the tangential 
component of the electric field is continuous, resulting in 

(V X B)i ^ (VxB)2 _ ^3g^ 

(Ti (72 



It follows from (g^) that the currents are much stronger 
in the inner layer in a quasi-steady state. Another way to 
see that the currents are pumped down the crustal area, 
is directly from (|l^): — —dxT] ^ d^cr / cP' ■ As the con- 
ductivity increases inwards, this part of velocity results in 
pushing down the magnetic flux. Therefore, if the initial 
currents are evenly distributed in the crustal area, the 
upper currents dissipate in short time scale (30 years), 
currents in deeper layers dissipate over longer time scales 
3000 years), while the whole crustal field lasts for a 
few million years. 

As we mentioned before, it is not known which part of 
the currents supporting the poloidal field is situated in 
the crust In any event, that part of the crustal 

currents can dissipate via our mechanism in a very short 
time scale while the field anchored in the core may remain 
for time scale comparable with the age of the universe. 
It is possible that pulsars with relatively low observed 
magnetic fields indicate a core component of ~ 10® G, 
while pulsars with fields of order 10^^ G are younger and 
have not had time to lose their crustal field component. 
As isolated pulsars lose their crustal magnetic field due 
to rapid decay, they also slow down, in the process cross- 
ing the death line to become unobservable. We suggest 
that as neutron stars in binary systems lose their crustal 
magnetic fields, they permit an increased rate of accre- 
tion that spins them up to give rise to the millisecond 
pulsar population. 
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